void TOTERG(Ener, Vir)
{
	double xi, yi, zi, Ener, eni, CORU, viri, Vir, rho;
	int i, jb;

	Ener = 0;
	Vir = 0;
	for (i = 1; i < NPART; i++)
	{
		xi = X(i);
		yi = Y(i);
		zi = Z(i);
		jb = i + 1;
		ENERI(xi, yi, zi, i, jb, eni, viri);
		Ener = Ener + eni;
		Vir = Vir + viri;
	}
	if (TAILCO) {
		rho = NPART/(BOX*BOX*BOX);
		Ener = Ener + NPART*CORU(RC, rho);
	}
}
